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CRITICALITY EXPERIMENTS AND ANALYSIS OF MOLYBDENUM -REFLECTED 
CYLINDRICAL URANYL -FLUORIDE - WATER SOLUTION REACTORS 


by Daniel Fieno, Thomas Fox, and Robert Mueller 
Lewis Research Center 

SUMMARY 

Clean criticality data were obtained from molybdenum -reflected cylindrical uranyl- 
fluoride - water solution reactors. Using ENDF/B molybdenum (Mo) cross sections, a 
nine energy group two-dimensional cylindrical (r,z) S n calculation of a reflected reactor 
configuration having a hydrogen to uranium -235 atom ratio of 319.4 predicted criticality 
to within 7 cents of the experimental value. One -dimensional cylindrical radial S n cal- 
culations were not adequate for predicting the criticality of these reactors; for the seven 
reactor configurations studied, the calculations were from 0.45 to 1.84 dollars too reac- 
tive. 

For these reactors, it was necessary to compute the reflector effective resonance 

integral by a detailed transport calculation at the core -reflector interface volume in the 

energy region of the two dominant resonances of natural Mo. For the 44.7-eV resonance 
95 

of Mo , the resonance integral was computed to be 7.20 barns while for the 131.3-eV 
96 

resonance of Mo , the resonance integral was computed to be 3.55 barns. 

The calculated fundamental prompt-mode decay constants at delayed critical were in 
good agreement with the values measured for these reactors. 

Molybdenum cross sections from the ENDF/B number 1111 evaluation are adequate 
for predicting the multiplication factor of the seven reactors presented in this study. 


INTRODUCTION 

The conceptual design of a compact fast reactor for space power application and a 
technology program to support the reactor are in progress at Lewis Research Center 
(ref. 1). The reactor uses uranium nitride as fuel that is clad and supported by a tanta- 

rt 

lum (Ta) alloy T-lll; lithium-7 (Li ) is used as coolant and the core is reflected by the 
molybdenum (Mo) alloy TZM, both radially and axially. The nuclear design and criti- 



cality experiments for the reactor are described in reference 2. 

It has been our experience that attempts to calculate the neutron multiplication fac- 
tor of past fast critical assemblies fueled with oralloy and containing or reflected by the 
refractory materials Ta and Mo, have yielded values significantly larger than the ex- 
perimental factor of unity (ref. 3). The bulk of these discrepancies are attributed to 
uncertainties in fast neutron cross sections; smaller portions of the discrepancies have 
been assigned to the use of fewer neutron groups and lower order two-dimensional 
transport analysis. 

In order to test the validity of available cross section evaluations, a series of inte- 
gral experiments have been performed employing spheres of refractory materials with 
central neutron sources (ref. 4). Measured and calculated sphere leakage spectra in the 
energy range 0. 5 to 11 MeV were compared as a test of the evaluated total, inelastic 
scattering, and reaction cross sections. Data for Mo showed general agreement but sig- 
nificant differences were observed for spheres of Ta, tungsten (W), and beryllium (Be). 
In order to obtain additional data and to widen the neutron energy range of comparison, a 
series of criticality experiments have been obtained using the NASA enriched uranyl- 
fluoride - water (UO 2 F 2 -H 2 O) cylindrical solution reactor (ZPR-I) in a configuration 
having a thick molybdenum radial reflector . Both criticality data and radial neutron 
leakage spectra for this reactor have been obtained. The present report presents an 
analysis of the criticality data. 

Because of the homogeneity of the cylindrical zones, the critical experiments pre- 
sent relatively clean configurations for analysis by both one- and two-dimensional, mul- 
tigroup, discrete -ordinate transport theory computer codes. The adequacy of the most 
recent ENDF/B Mo cross section compilation may be tested by these calculations. 


REACTOR AND EXPERIMENT PROCEDURES 

The ZPR-I UOgFg-HgO solution reactor used in these tests is a 25. 4 -centimeter - 
diameter cylindrical core tank that is reflected radially with a 10. 14 -centimeter -thick 
Mo annular reflector . The experimental program involved the measurement of two basic 
parameters; namely, critical solution height and the fundamental prompt-mode decay 
constant. Brief descriptions of the critical experiments and the pulsed source experi- 
ments are presented. More complete descriptions of the criticality and pulsed source 
experiment procedures are presented in references 5 and 6, respectively. 


ZPR-I with Radial Molybdenum Reflector 

ZPR-I is a solution reactor that uses an aqueous solution of highly enriched (93.2 
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percent uranium-235 (U 0 )) uranyl fluoride (UOgFg) as the fuel that is readily varied in 
concentration. The fuel solution is stored in geometrically safe storage tanks and the 
critical configuration is achieved by remotely transferring the fuel solution from storage 
to the stainless-steel reactor tank. In all cases where a critical configuration was to be 
achieved, the reflector arrangement was assembled first and then the fuel concentration 
was adjusted so that the core and reflector heights were the same at delayed critical. 

The physical arrangement of the core and reflector is shown in figure 1. The radial 
dimensions are presented in table I. The reactor configuration is located at the approxi- 
mate center of a large room and is elevated about 20 centimeters above a working grid 
platform to minimize interactions with nearby structure. The Mo reflector blocks were 
also stacked on a cylindrical spool -like aluminum platform designed to minimize the 
structural interactions. 

The Mo blocks were fabricated by pressing and sintering and were machine finished. 
The material is 99.95 percent Mo and has an average density of 9.88 grams per cubic 
centimeter (96. 8 percent of theoretical density). The reflector annulus has been divided 
into eight equal segments and is constructed from as many as 32 blocks which are 20. 32 
centimeters in height and eight blocks which are 9. 87 centimeters in height. The size of 
the blocks permitted manual handling and reflector heights could be constructed to a 
maximum of 91 centimeters in approximately 10-centimeter increments. The radial 
thickness of the blocks in the radial dimension is 10. 14 centimeters and when stacked 
in a tight array form a cylindrical chimney with an inside diameter of 26. 87 centimeters 
and an outside diameter of 47.15 centimeters. 

The reactor tank is constructed of type 304 stainless steel, 0.16 centimeter thick, 
and has an outside diameter of 25.72 centimeters. Allowing for the corrosion protective 
coating on the tank wall, there is a gap of 0. 57 centimeter between reactor tank and re- 
flector when only the blocks are used as a reflector. This gap was designed to accom- 
modate the insertion of a cadmium (Cd) sleeve between the core and reflector as shown 
in figure 1. The Cd has been constructed of 0.0 89 -centimeter -thick Cd sheet canned in 
stainless steel. The inserts are in the form of half cylinders constructed to eliminate 
gaps in the Cd layer that would permit radial leakage of subcadmium neutrons. The Cd 
sleeve was adjusted in length to match the height of the reflector used. Molybdenum in- 
serts about 0.48 centimeter thick were also fabricated to fill the radial gap. Owing to 
dimensional tolerances, this is about the maximum amount of material that can be 
readily fitted into the gap and fills about 85 percent of the gap. These pieces were also 
fabricated in half cylinders of enough different sizes to permit a match of all of the re- 
flector heights used during this study. The density of this material was measured and 
found to be within 1 percent of the theoretical density of 10.2 grams per cubic centi- 
meter . 
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Figure 1. - Schematic diagram of NASA ZPR-I solution reactor with annu- 
lar molybdenum reflector. (Cadmium insert sleeve is shown between 
stainless-steel tank wall and molybdenum reflector. ) 


TABLE I. - RADIAL DIMENSIONS OF THE NASA 
ZPR-I SOLUTION REACTOR HAVING AN 
ANNULAR MOLYBDENUM 
REFLECTOR 



Inside radius, 
cm 

Outside radius, 
cm 

Fuel solution a 

0 

12.70 

Type 304 stainless - 

12.70 

12.86 

steel tank 



Insert region b 

12.86 

13.43 

Mo reflector 

13.43 

23.57 


a U 0 2 F 2 salt in HgO solution with concentration varied 
so that fuel solution height was same as Mo reflec- 
tor height for given experiment. 

^Insert region either was left void or contained 0.089- 
cm-thick Cd sleeve clad in stainless steel, or 0.48- 
cm Mo sleeve. 


Pulsed Neutron Source Experiments 


The experimental arrangement used to conduct the pulsed neutron source experi- 
ments is the same as described in reference 6. The pulsing head from the neutron gen- 
erator is located on the axis of the reactor core directly beneath the tank and the neutron 
detector is located in the core at a position found to be optimum for studying the funda- 
mental prompt mode (see fig. 1). The neutron pulse generator is a small positive-ion 
sealed tube accelerator which produces 14.3-MeV neutrons by means of the ^H(d,n)^He 
reaction. The neutron pulse intensity is about 10 neutrons per pulse with the pulse 
width being 3. 5 microseconds and the rate about 5 pulses per second. The neutron de- 

Q 

tector probe is a small 1/v gas proportional counter filled with helium-3 ( He) gas. 
This probe is less than 1 centimeter in diameter and has a sensitive volume of 1 cubic 
centimeter which approximates a point detector. The pulses from the detector were fed 
through a preamplifier and amplifier to a 400-channel analyzer which was operated in a 
multiscaler mode using 100- and 200 -microsecond channel widths. System resolution 
time has been measured as 1.12 microseconds with this arrangement. 

The output data from the 400 -channel analyzer is reduced using the digital computer 
program GRIPE n developed by Kaufman (ref. 7). This program calculates the funda- 
mental prompt-mode decay constant and can be used to calculate the shutdown reactivity 
of a subcritical solution system. The data reduction methods are described in detail in 
reference 6. 


Experimental Program 

A series of experiments were made using the Mo reflected solution reactor. Sev- 
eral measurements were made of the reactivity effect associated with the voided gap be- 
tween the core and the reflector, with the Cd sleeve inserted, and with the Mo sleeve in- 
serted. Fuel concentration was varied over a range of hydrogen to uranium-235 atom 
ratios H/X of 319.4 to 615 during the tests. Delayed critical configurations were 
achieved with nominal reflector heights of 30. 2, 40.6, 50.5, and 70. 8 centimeters. The 
maximum height was limited by the height of the core tank to 75 centimeters. 

Two sets of delayed critical experiments were run at each reflector height. The 
first set of experiments had the 0. 57-centimeter gap between the core and reflector 
voided; for the second set the gap was filled with Mo inserts. To achieve the same crit- 
ical height for the two sets required an adjustment in fuel concentration. The worth of 
the Mo sleeve can thus be measured in terms of a change in fuel concentration. If the 
reactivity is known as a function of fuel concentration, it is thus possible to get a reac- 
tivity worth for the Mo sleeve by this technique. 


E-6812 
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Another measurement of the reactivity worth of the Mo insert as well as the worth 
of a Cd insert was made using the pulsed neutron source technique. The most reactive 
configuration, which was the arrangement with the Mo insert, was used as the reference 
case. By using the fuel concentrations for the clean core cases previously run, the de- 
layed critical condition was established with the pulsing equipment in place for fuel con- 
centrations with H/X values of 330, 540, and 615. At each fuel concentration there 
was a small increase in critical height due to the presence of the detector probe and the 
proximity of the pulsing head. The new critical height was used throughout the pulsed 
experiments at each fuel concentration. The Cd insert and void gap arrangements were 
both subcritical at the reference fuel heights for all fuel concentrations used . 

The reactivity worth measured by this pulsed source technique is considered to be 
a better value than the worth that could be obtained as a function of fuel concentration 
change (see ref. 6) . Most important in this application, however, is the fact that it is 
not necessary to measure the fuel concentration to obtain the reactivity worth. It is 
only necessary to assure that the fuel concentration has not changed between the critical 
case and the shutdown case being compared. Although fuel concentration is a difficult 
variable to measure with precision, it is relatively easy to assure that no significant 
change in concentration has occurred in repeating a reference case. The pulsed source 
technique used in this way is most advantageous. 


TRANSPORT CALCULATIONS 

Calculation of the Absolute Criticality of the 
Molybdenum Reflected Solution Reactors 

The reactor configurations considered are such that a two-dimensional cylindrical 
(r, z) calculation accurately describes the geometry. This means that the absolute mul- 
tiplication factor of each reactor configuration could be obtained subject only to errors 
in the neutron cross section data if parameters such as the number of neutron energy 
groups, which are controllable, are adequate. Thus, a two-dimensional calculation for 
a given reactor configuration could serve as a standard for assessing the error in the 
critical height obtained from a corresponding one -dimensional radial calculation. 

Since the delayed critical reactors were achieved with essentially the same fuel so- 
lution heights and Mo reflector heights, two-dimensional effects are minimized. This 
implies that one -dimensional cylindrical radial calculations with an appropriate axial 
buckling could be done for the large number of cases required and that the absolute mul- 
tiplication factor of a given reactor could be obtained by applying a correction factor 
based on the standard two-dimensional result. However, reactor parameters, for ex- 
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ample, the effective delayed neutron fraction which are needed in the analysis of the 
pulsed source experiments are not sensitive to small errors in the absolute reactivity 
of a given reactor . 

A two-dimensional cylindrical (r,z) S n calculation using the DOT code (ref. 8) was 
made to obtain the absolute multiplication factor for one of the delayed critical reactors. 
This reactor had a hydrogen to uranium-235 atom ratio of 319. 4. The experimental re- 
actor (and reflector) height was 30.74 centimeters. This reactor was critical with the 
gap between the Mo reflector and the stainless -steel core vessel left void. 


One-Dimensional S n Calculations 

The bulk of the calculations for this work were performed using a one -dimensional 
S n transport code developed at the Lewis Research Center. Radial calculations were 
performed for the cylindrical critical systems with the fuel solution, stainless -steel con- 
tainment vessel, gap (voided or filled with Mo or Cd) between the containment vessel and 
the reflector, and the Mo reflector explicitly included in the calculations. The leakage 
through the ends of the critical cylinders was accounted for through a buckling term. 

The eigenvalue for the system was taken to be the height of these cylindrical reactors 

P 

for a given value of the fundamental prompt-mode decay constant q!q . 

The calculational procedure followed is described in reference 6 and is summarized 

P 

briefly here. One -dimensional, radial S n calculations with an assigned value of aig 
were made for a given critical experiment to obtain both the fundamental prompt-mode 
flux and the static -mode adjoint flux. These calculations yield the reactor height for the 
assigned value of and the calculated reactivity Pg($) by using either of the follow- 
ing two equivalent expressions: 


„ K , f - 1 
c / a>\ ett 


P*($) = 


K 


eff ^0 


or 


p£($)=4r- +1 

fo 

An 


c 

where p g ($) is the calculated value of the static reactivity in dollar units, is the 


static -mode multiplication factor, /3 q is the fundamental prompt-mode effective delayed 

neutron fraction, A n is the fundamental prompt-mode generation time in seconds, and 
P u 

a?Q is the fundamental prompt -mode decay constant in inverse seconds. The value of 

a? is varied (and correspondingly , the reactor height) until the reactivity p^($) is 
vr c S 

zero. This calculated height for p s ($) equal to zero gives the calculated delayed criti- 
cal configuration which is compared to an experimental value. 

Nine neutron energy groups were used, eight fast and one thermal group which in- 
cluded an up-scattering transfer component. The fast group neutron cross sections for 
the UOgFg-HgO fuel solutions, the type 304 stainless-steel reactor tank, and Cd were 
obtained using the GAM n code (ref. 9). The neutron cross sections for the Mo reflec- 
tor were obtained from the ENDF/B data files and averaged over the fundamental mode 
fluxes and currents generated by GAM n for the various fuel solutions. The thermal 
group cross sections for all materials were obtained using the GATHER II code (ref. 10). 

The energy boundaries for each neutron energy group are given in table n. The delayed 
235 

neutron data for U which are needed for the calculations are given in table HI 

(ref. 11). All of the delayed neutrons are assumed to appear in energy group four where 

the energy varies from 0.22 to 0. 82 MeV. 

The one -dimensional radial S n transport calculations were all performed using an 
angular quadrature order with elastic scattering treated through the P.^ order. The 
work of reference 12, although performed for an unreflected solution reactor, indicates 
that increasing the angular quadrature order, the elastic scattering order, or the num- 
ber of thermal groups had a negligible effect on the eigenvalue determined for the de- 
layed critical data. 

TABLE IB. - DELAYED NEUTRON DATA FOR 


TABLE II. - NEUTRON ENERGY GROUP 


BOUNDARIES 


Group 

Energy 

Lethargy 

1 

14.92 to 4.966 MeV 

-0.4 to 0.7 

2 

4.966 to 2.231 MeV 

0.7 to 1.5 

3 

2.231 to 0. 821 MeV 

1 . 5 to 2. 5 

4 

0.821 to 0.224 MeV 

2.5 to 3.8 

5 

0.224 MeV to 9.12 keV 

3.8 to 7.0 

6 

9.12 keV to 454 eV 

7.0 to 10.0 

7 

454 to 8.32 eV 

10.0 to 14.0 

8 

8.32 to 0.414 eV 

14.0 to 17.0 

9 

0.414 to 0 eV 

17.0 to 00 


URANIUM-235 


[From ref. 11, p. 90; /3 = 0.0065.] 


Delayed 

group 

Decay 

constant, 

sec™ 1 

Relative 

abundance, 

Bj/0 

Delayed neutron 
fraction, 

B i 

1 

0.0124 

0.033 

2.145X10 -4 

2 

.0305 

.219 

14.235X10 -4 

3 

.111 

.196 

12.740X10 -4 

4 

.301 

.395 

25.675X10 -4 

5 

1.14 

.115 

7.475X10" 4 

6 

3.01 

.042 

2.730X10" 4 
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Calculated kinetic parameters are presented in table IV for four fuel solutions hav- 
ing H/X values of 319.4, 447, 528, and 596. These fuel concentrations, whose atom 
densities are given in table V, were chosen so that fuel solution height and the Mo side 
reflector height was as nearly the same as possible for four critical experiments. The 
nominal experimental heights for these cylindrical systems were 30. 2, 40. 6, 50. 5, 
and 70. 8 centimeters. These calculations were performed for the gap between the core 
tank and Mo annular side reflector left void. 

Figure 2 shows the variation in the fundamental prompt-mode effective delayed neu- 
tron fraction /3_ with respect to H/X for the fuel solution. The value of 8 , the delayed 
^ 235 

neutron fraction, was taken to be 0.0065 for U . Figure 3 shows the fundamental 

prompt-mode generation time Aq as a function of H/X. Figure 4 shows the ratio of 

the parameter /3 q/A q as a function of H/X. Since cyq = -/3 q/Aq at delayed critical, 

the negative of the data shown in figure 4 is also the fundamental prompt-mode decay 

constant q!q at delayed critical. Figure 5 shows the reactivity coefficient Ap^($)/Ah 

as a function of H/X. This coefficient gives the worth in dollars per centimeter change 

in reactor height at delayed critical. Figure 6 shows the reactivity coefficient Ap C ($)/ 
p s 

Aq?q as a function of H/X. This coefficient gives the change in reactivity due to a 
change in the fundamental prompt-mode decay constant. 


TABLE IV. - TABULATION OF CALCULATED DELAYED CRITICAL DATA 

[Calculated values are for void gap between stainless-steel core tank and molybdenum 
annular reflector. | 


Hydrogen to uranium -235 atom ratio, 

H/X 

319.4 

447 

528 

596 

Effective delayed neutron fraction, 
^0 

7.907X10' 3 

7.752X10 -3 

7.652X10 -3 

7.564X10' 3 

Neutron generation time, 
A 0 , ju-sec 

26.84 

35.27 

40. 15 

44.00 

Ratio of effective delayed neutron frac- 
tion to neutron generation time, 

^q/A 0’ sec ~* 

294.6 

219. 8 

190.6 

171.9 

Reactivity coefficient, 
AK eff /Ah ’ cm_1 

9. 896X10 " 3 

5. 207X10 -3 

2.923X10' 3 

1. 408X10" 3 

Reactivity coefficient, 
A Pg($)/Ah, {-cm' 1 

1.252 

0.672 

0.382 

0. 186 

Reactivity coefficient, 
Apg( $)/A c*q , $-sec 

3. 394X10 -3 

4. 550X10" 3 

5. 247X10' 3 

5. 817X10" 3 
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TABLE V. - ATOM DENSITIES OF FUEL SOLUTIONS REQUIRED TO 


YIELD NOMINAL EXPERIMENTAL DELAYED CRITICAL HEIGHTS 
OF 30.2, 40.6, 50.5, AND 70 . 8 CENTIMETERS 
[Void gap between stainless -steel tank and molybdenum reflector.] 


Nuclide 

Hydrogen to uranium 

-235 atom ratio, H/X 


319.4 

447 

528 

596 


Delayed critical height, h, cm 


30.2 

40.6 

50.5 

70.8 


Atom densities, atoms/b-cm 

Hydrogen 

0.066019 

0.066233 

0.066315 

0.066367 

Oxygen 

.033453 

.033435 

.033427 

.033422 

Fluorine 

4.436X10' 4 

3.179X10' 4 

2.695X10" 4 

2. 3 89X10 “ 4 

Uranium-23 5 

2.067X10' 4 

1.482X10' 4 

1.256X10 -4 

1.114X10 -4 

Uranium-238 

1. 506X10' 5 

1.079X10' 5 

CD 

cn 

X 

H-*- 

o 

1 

Oi 

8.1 1x10 " 6 


8. Ox 10’ 3 



UJ 


7. 5l L I I [ I I 

300 350 400 450 500 550 600 

Hydrogen to uranium -235 atom ratio, H/X 

Figure 2. - Effective delayed neutron fraction at delayed critical as a 
function of hydrogen to uranium-235 atom ratio for NASA solution 
reactor with molybdenum side reflector. 
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Hydrogen to uranium-235 atom ratio, H/X 


Figure 3. - Neutron generation time at delayed critical as a function 
of hydrogen to uranium -235 atom ratio for NASA solution reactor 
with molybdenum side reflector. 



Hydrogen to uranium-235 atom ratio, H/X 


Figure 4. - Ratio of effective delayed neutron fraction to neutron gen- 
eration time at delayed critical as a function of hydrogen to uranium- 
235 atom ratio for NASA solution reactor with molybdenum side re- 
flector. 
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Reactivity coefficient, Ap c ($)/Aaj?, $-sec 


Hydrogen to uranium-235 atom ratio, H/X 

Figure 5. - Reactivity coefficient in dollars per centimeter of reactor 
height at delayed critical as a function of hydrogen to uranium- 
235 atom ratio for NASA solution reactor with molybdenum side 
reflector. 


6. OxlO' 3 



Hydrogen to uranium -235 atom ratio, H/X 


Figure 6. - Reactivity coefficient in dollar-seconds at delayed critical 
as a function of hydrogen to uranium-235 atom ratio for NASA so- 
lution reactor with molybdenum side reflector. 





MOLYBDENUM CROSS SECTIONS 


ENDF/B Molybdenum Cross Section Evaluation 

The ENDF/B Mo data was compiled in 1966 by Pennington and Gajniak (ref. 13) and 
issued as ENDF/B material number 1025. The primary source of information for this 
evaluation by Pennington and Gajniak is the evaluation prepared by Schmidt (ref. 14). 
Pennington and Gajniak made the following modifications to the data compiled by Schmidt: 

(1) Resonance parameters were taken from reference 15, plus some smooth capture 
cross sections were used from 4 eV to 1 keV. 

(2) Inelastic scattering cross sections from 200 keV to 1. 5 MeV, for four resolved 
levels, were obtained from reference 16. 

(3) The (n,2n) and (n,3n) cross sections were obtained from reference 17. 

(4) The data range was extended from 10 to 15 MeV: (1) the total cross section by 
using the data in references 18 and 19, (2) the (n,y) cross section by assuming a 1/E 
behavior, and (3) the nonelastic cross section was assumed to be constant. 

(5) Above the (n,2n) threshold, the inelastic cross section was assumed to be the 

difference between the nonelastic cross section and the sum of the cross sections for 

\ 

various other reactions. 

(6) Nuclear temperatures were estimated for continuum inelastic scattering as in 
reference 20. 

(7) Legendre coefficients for elastic scattering were obtained up to 6 MeV from ref- 
erence 21 and above 6 MeV, the coefficients were obtained from optical model calcula- 
tions. 

The Mo cross sections have been further revised and issued as ENDF/B material 
number 1111. A number of changes in the ENDF/B #1025 evaluation were made. The 
most important change for the purpose of this study was in the high energy (n,y) cross 
sections. This (n,y) cross section data was reduced in value so that computed central 
reactivity worths in fast assemblies correspond to experimental reactivity determina- 
tions . 


Result of an Integral Experiment to Test Inelastic Scattering 
Cross Sections of Molybdenum 

In order to test the validity of the ENDF/B #1111 Mo cross sections, an integral 
experiment has been performed employing a Mo sphere with a centrally located neutron 
source (ref. 4). Here, measured and calculated sphere leakage spectra in the energy 
range 0. 5 to 11 MeV are compared as a test of the evaluated inelastic scattering cross 
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sections. A 54 Ci Am-Be source was used as a source of fast neutrons at the center of 
a 22. 96- centimeter-diameter Mo sphere whose shell thickness was 8. 38 centimeters. 

The Mo was in the form of a powder having an average density of 3. 73 grams per cubic 
centimeter. This gave a shell thickness of 0. 74 total mean free paths at 4 MeV. The 
leakage spectrum was measured at a distance of 200 centimeters using an NE213 liquid 
scintillator neutron spectrometer. Room scattered neutron backgrounds were accounted 
for by using a shadow cone. The leakage spectrum was calculated with 51 energy groups • 

and an S n discrete ordinate transport code. 

The results of this integral experiment are shown in figure 7. The agreement be- , 

tween experiment and calculation is satisfactory and similar to the results obtained for 
the bare source. From this integral experiment, one can conclude that the inelastic 
scattering and total cross sections for Mo in the energy range from 0. 5 to 11 MeV are 



Figure 7. - Leakage spectrum for source enclosed by molybdenum shell (54 Ci Am-Be neutron source with intensity of 1. 3x10® neutrons/sec). 
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satisfactory. However, the other Mo ENDF/B # 1111 neutron cross sections need to be 
further tested by experiments such as the criticality experiments presented in this work. 


Resonance Cross Sections 

Molybdenum has a number of absorption resonances which yield an infinitely dilute 
resonance integral of 29. 3 barns for the ENDF/B #1111 evaluation. Table VI lists the 
resonances which are the largest contributors to the resonance integral. The 44. 7-eV 
resonance of Mo^ and the 131. 3-eV resonance of Mo^ contribute 22. 8 barns to the in- 
finitely dilute resonance integral. The remaining resonances of Mo give small resonance 
integrals and it is expected that the energy and spatial self -shielding factors for these 
resonances will not differ appreciably from unity. 

Since criticality calculations using Mo absorption cross sections computed using the 
GAM II approximate resonance calculation gave reactivity differences of $3. 76 in the 
critical heights as compared to absorption cross sections based on the infinitely dilute 
resonance integral, an accurate evaluation of the resonance integral of the two main 
resonances of Mo was needed. Therefore, effective resonance integrals for the two 
most important resonances of Mo were computed using the following procedure: For 
each resonance considered, Doppler broadened absorption and scattering cross sections 


TABLE VI. - RESONANCES WHICH ARE THE MOST 
IMPORTANT CONTRIBUTORS TO THE INFINITELY 


DILUTE MOLYBDENUM RESONANCE INTEGRAL 3 


Isotope 

Resonance energy, 
eV 

Resonance integral 
contribution for 
natural molybdenum, 
b 

Mo 95 

44. 7 

17.3 


358. 5 

.4 

Mo 96 

131.3 

5. 5 

Mo 97 

70.9 

0.5 

Mo 98 

12.1 (p wave) 

0. 8 


467.4 

.7 

Mo 100 

363.7 

0.6 


a Infinitely dilute resonance integral for ENDF/B 
1111 Mo data is 29.3 b. 
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were computed for a temperature of 293 K. Then down-scattering transfer cross sec- 
tions Pq and Pj were computed for Mo and water for the group structure chosen to 
represent each resonance. Finally, P^S^ transport calculations were performed for 
the actual reactor geometry with water representing the fuel solution. The results of 
these transport calculations for the effective resonance integral are presented in 
table VP. 

The data in table VII indicate that the calculation of the effective resonance integral 
was not sensitive to the geometry used or to the height of the configurations Calcula- 
tions were also made in periodic slab geometry for the same thickness of water but 
varying the Mo thickness. The effective resonance integral for both the 44.7- and 
131.3-eV resonances approached the infinitely dilute values as the Mo thickness became 
progressively smaller. 

For the Mo reflected solution reactors considered in this report, the effective reso- 

95 

nance integral of the 44.7-eV resonance of Mo was taken as 7.20 barns, while the ef- 
fective integral of the 131. 3-eV resonance of Mo was taken as 3.55 barns. Both reso- 
nances occur in group 7 of the energy group structure presented in table II. The 
ENDF/B # 1111 Mo absorption cross sections in group 7 were adjusted so that the effec- 
tive resonance integral of the two resonances considered was 10.75 barns. 


TABLE Vn. - EFFECTIVE RESONANCE INTEGRALS FOR THE 44.7-eV 
RESONANCE OF MOLYBDENUM-95 AND THE 131.3-eV 
RESONANCE OF MOLYBDENUM-96 


[The calculations were using a water region 12. 7 cm thick and a 

molybdenum thickness of 10. 14 cm. ] 


Resonance energy, eV 

44.7 (Mo 95 ) 

131.3 (Mo 96 ) 

Inifinitely dilute resonance integral, b 

17.3 

5.5 

Number of energy groups in S 

26 

21 

calculation 



Computed resonance integral, b 



Cylindrical geometry, 30.45-cm height, 

7.20 

3.55 

no return current boundary condition 



Slab geometry, 30.45-cm height, no 

7. 15 

3.54 

return current boundary condition 



Periodic slab geometry, 30.45-cm 

7.14 

3.54 

height 



Periodic slab geometry , infinite height 

7. 16 
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COMPARISON OF EXPERIMENTS AND CALCULATIONS 


The comparison of the calculations with experiment will be made in terms of two 
quantities: (1) the delayed critical height, and (2) the fundamental prompt-mode decay 
constant. These two quantities were measured for each reactor configuration and are 
sufficient to establish the absolute criticality. 


Delayed Critical Heights 

Void gap between stainless-steel core tank and molyb denum reflector . The de- 
layed critical height was measured for four fuel solution concentrations. These fuel so- 
lutions had hydrogen to uranium 235 atom ratios H/X of 319.4, 447, 528, and 596 which 
gave critical reactor heights of 30.74, 40.54, 50.55, and 70.42 centimeters, respec- 
tively. The experimental heights were corrected to a fuel solution temperature of 20° C 
and for loss of water by evaporation. Both of these corrections were small in compari- 
son with the reactor heights . 

The reactor with the smallest delayed critical height, that is, the reactor with a fuel 
solution having an H/X of 319.4, was computed using the two-dimensional transport 
program DOT. This calculation served as a standard for assessing the error in the ab- 
solute criticality determination obtained from a corresponding one-dimensional S n 
transport calculation. This DOT calculation was made using cylindrical (r,z) geometry 
and used nine energy groups, a scattering order, and an S 4 angular quadrature . 

The experimentally determined height of 30.74 centimeters was used in this calculation. 
The computed effective static multiplication factor was 1.00054 which indicates, upon 
using the value of /3 eff for this H/X from table IV, that the calculation is within 
6. 8 cents of the experimental multiplication factor of unity. 

A one -dimensional cylindrical radial calculation for this same reactor yielded a de- 
layed critical height of 29.27 centimeters. Using the reactivity coefficient Ap c ($)/Ah 
for this H/X from table IV indicates that the one -dimensional calculation is too reac- 
tive by 1.84 dollars. Similar results can be expected for the other three critical 
reactors. 

A comparison of experimental and one -dimensional calculated values of the delayed 
critical height, h, as a function of H/X for the four fuel solutions considered has been 
made in figure 8 and table VIII. All of these calculations predicted critical heights which 
were less than the experimental delayed critical heights. As shown in table VIII, these 
one -dimensional calculations were over-reactive by amounts of 184 to 115.4 cents for 
solution concentrations which varied from an H/X of 319. 4 to 596. 
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Figure 8. - Delayed critical height as a function of hydrogen to 
uranium-235 atom ratio for NASA solution reactorwith molybde- 
num side reflector. (Void gap between stainless-steel tank and 
molybdenum reflector.) 


TABLE VIE. - TABULATION OF DELAYED CRITICAL HEIGHTS AND 
FUNDAMENTAL PROMPT-MODE DECAY CONSTANTS 


[Computed values are based on one -dimensional cylindrical radial 
calculations. ] 


Hydrogen to 
uranium -2 3 5 
atom ratio, 
H/X 

Delayed critical height, 
h, 

cm 

Fundamental prompt -mode 

decay constant, 

P 

a 0’ 

sec"* - 

Experimental 3 

Calculated 

Experimental 

Calculated 


Void gap between stainless- 

-steel tank 


319.4 

30. 7±0. 1 

29.27 ( b 184. 0) 

-281±2 

-293 

447 

40. 5±0. 1 

38.02 ( b 169 . 3) 

-214±2 

-218 

528 

50. 6±0. 1 

47. 53 ( b 1 1 5 . 4) 

-187±2 

-191 

596 

70. 4±0. 1 

61.97 ( b 157. 2) 

-166±2 

-173 


Gap filled with molybdenum 


330 

30. 2±0. 1 

29.6 ( b 66. 9) 

-274±2 

-285 

540 

49. 9±0. 1 

47. 8 ( b 71. 8) 

-185±2 

-186 

615 

68. 9±0. 1 

65.6 ( b 45. 8) 

-165±2 

-166 


Experimental heights corrected to a fuel solution temperature of 20° C 
and for evaporation. 

^Reactivity worth in cents of height difference between experimental and 
calculated values. 


Molybdenum sleeve between stainless-steel tank and molybdenu m reflector . - Three 
critical reactor experiments were performed for a Mo sleeve filling the 0. 57-centimeter 
gap between the stainless-steel reactor tank and the Mo reflector. A comparison of ex- 
perimental and one -dimensional calculated values of the delayed critical height, h, as a 
function of H/X is made in figure 9. These data are also presented in table VIII. 

These critical experiments were performed for slightly leaner fuel concentrations than 
those with the voided gap. The leaner fuel concentrations were necessary to maintain 
the fuel solution height the same as the Mo reflector height for these more reactive ex- 
periments . 



Figure 9. - Delayed critical height as a function of hydrogen to uranium-235 
atom ratio for NASA solution reactor with molybdenum side reflector. (Gap 
between stainless-steel tank and molybdenum reflector filled with molybde- 
num.) 


At an H/X of 330, the calculated height is 0. 56 centimeter less than the experi- 
mental height. At H/X of 540 and 615, the calculated heights are 2.08 and 3.27 centi- 
meters less than the measured height. By using the calculated reactivity coefficient 
Ap^($) /Ah given in table IV, the difference of the calculated and experimental delayed 
critical heights in terms of reactivity is about 67, 72, and 46 cents for the reactors with 
H/X of 330, 540, and 615, respectively. 


Fundamental Prompt -Mode Decay Constants 

p 

The fundamental prompt-mode decay constant (Xq has been measured for the vari- 
ous critical systems considered by using the pulsed source technique. The measured 
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data are shown in figure 10 as well as being presented in table VHI. The decay constant 

p 

o/q is primarily, for these reactors, a function of the fuel concentration. This can be 
seen to some extent in figure 10 since the experimental data all fell on one curve regard- 
less of whether or not the gap between the core tank and reflector is left void or filled 
with a Mo sleeve. 

p 

A calculated curve is also shown in figure 10. The calculated values of o!q are in 
good agreement with the measured values. This can be best judged by using the reac- 
tivity coefficient Ap^( $) /&(Xq as given in table IV, which indicates that differences in 


Experimental data (void yap between stainless- 



Figure 10. - Fundamental prompt-mode decay constant at delayed critical as a 
function of hydrogen to uranium-235 atom ratio for NASA solution reactor 
with molybdenum side reflector. 


P 

(Jq correspond to differences in reactivity of less than 5 cents over the range of fuel 
concentrations used. Physically, this means that the nuclear cross section data used 
provide a good representation of the transport, moderating, and thermalization proper- 
ties for the fuel solution. 


Reactivity Worth of Molybdenum Sleeve and Cadmium Insert 

Pulsed source subcritical measurements provide an excellent means for determining 
the reactivity worths caused by removal of the Mo liner or its replacement by a 
stainless -steel clad Cd sleeve for the three experiments for which the reactors were 
critical with the Mo filled gap. If the experiments are done in this manner, the fuel so- 
lution height is the same as the Mo reflector height and thus, two-dimensional chimney 

P 

effects are avoided. Experimental measurements of q/q are presented in table IX for 
various configurations. 
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TABLE IX. - EXPERIMENTAL FUNDAMENTAL PROMPT-MODE DECAY 


CONSTANTS FOR MOLYBDENUM REFLECTED SOLUTION REACTORS 


Experimental configuration 

Hydrogen to uranium-235 atom ratio, 



H/X 



330 

540 

615 


Fundamental prompt-mode decay constant, 
Oq, sec~* 

Gap filled with molybdenum 

-274*2 

-185*2 

-165*2 

(delayed critical) 




Voided gap (subcritical) 

-604*4 

-380*4 

-320*2 


a (-l. 20) 

a (-l. 05) 

a (-0. 94) 

Stainless -steel clad cadmium 

-587*4 

-425*4 

-380*2 

sleeve (subcritical) 

a (-l. 14) 

a (-l. 30) 

a (-l. 30) 


a Subcritical reactivity worth in dollars relative to the reactor with gap 
filled with Mo. 


The method of Simmons and King (ref. 22) can be used to determine the experimen- 
tal subcritical reactivities for the configurations involved since the departure from the 
delayed critical configuration is not too great. By this method, the subcritical reac- 
tivity is given by 


($) = - 



+ 1 


where chq is the decay constant for the subcritical configuration and is the 

decay constant for the corresponding delayed critical configuration. Reference 5 gives a 
discussion of this method for determining subcritical reactivities. 

Removal of the Mo sleeve from the delayed critical reactors with fuel solution H/X 
of 330, 540, and 615 gives an experimental reactivity worth for the Mo sleeve of -1.20, 
-1.05, and -0.94 dollars, respectively. Calculated reactivity worth of the Mo sleeve 
was -39, -49, and -53 cents for the systems with an H/X of 330, 540, and 615, re- 
spectively. The calculated reactivities are lower than the experimental values by fac- 
tors of one-third to one-half. These data are shown in figure 11 along with the experi- 
mental data. 

Replacement of the Mo sleeve by a stainless -steel clad Cd liner gave experimental 
reactivity worths of -1.14, -1.30, and -1.30 dollars for the reactors with an H/X of 
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- O- Experimental reactivity worth with void between tank and reflector 
-O- Experimental reactivity worth with cadmium between tank and 
reflector 

Calculated reactivity worth of w id region between tank and 

reflector 



Figure 11. - Subcritical reactivity caused by removal of molybdenum sleeve 
as a function of hydrogen to uranium-235 atom ratio for NASA solution 
reactor with molybdenum side reflector. 


330, 540, and 615, respectively. Calculated reactivity worths as given in figure 11 
were -70, -89, and -97 cents for the corresponding values of H/X, respectively. The 
calculated reactivities are lower than the experimental reactivities by 39 to 25 percent. 

There are a number of reasons why the calculated reactivity is much lower than the 
experimental reactivity. The most important factor is the inability of the one- 
dimensional calculations to properly treat the neutron streaming in the voided gap. 

Here, the two-dimensional calculations would be more appropriate. Another factor is 
that the measured thickness of the various inserts and claddings varies within fabrica- 
tion tolerances and limitations. 


Thermal Flux Distributions 

The thermal neutron flux distribution in the radial direction at the reactor midplane 
was measured using dysprosium -aluminum foils for the reactor system having a fuel 
solution with an H/X of 319. 4. The dysprosium activity was measured for the insert 
region between the stainless -steel core tank and the Mo reflector voided and for an in- 
sert consisting of a Cd sleeve clad with stainless steel. These measured distributions 
are shown in figure 12 and are normalized to the activity measured at the center of the 
reactor . 
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Q Dysprosium foil activity for stainless-steel clad cadmium sleeve 
in insert region 

O Dysprosium foil activity for voided insert region 

Calculation for wided insert region 

Calculation for stainless-steel clad cadmium sleeve in insert 

region 



Figure 12. - Relative thermal flux distribution for reactor with hydrogen to ura- 
nium -235 atom ratio of 319.4. 


Calculated radial thermal flux distributions are also shown in figure 12. The mea- 
sured and calculated thermal flux distributions are in agreement with each other and ex- 
hibit the same features. The flux distributions for this reactor with and without the Cd 
insert are essentially the same within the core region. However, the flux distributions 
do differ, as expected, within the insert region. Thus the Mo reflector is effective pri- 
marily as a reflector of fast neutrons and has little effect on the thermal neutron flux 
distribution in the core because of its large thermal absorption cross section. 


CONCLUSIONS 

A number of conclusions result from this criticality study of molybdenum-reflected 
cylindrical uranyl -fluoride - water solution reactors: 

1 . A nine group two-dimensional cylindrical (r , z) transport calculation predicted 
the criticality to within 7 cents. 

2. Nine group one -dimensional cylindrical radial transport calculations for the 
seven reactor configurations studied were from 0.45 to 1.84 dollars too reactive. 

3. The effective resonance integral of molybdenum (Mo) must be accurately calcu- 
lated for these reactors. The effective resonance integral of the 44.7-eV resonance of 
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Mo was computed to be 7.20 barns, while the effective resonance integral of the 

98 

131. 3-eV resonance of Mo was computed to be 3. 55 barns. 

4. The calculated fundamental prompt -mode decay constants at delayed critical were 
in good agreement with the values measured for these reactors. 

5. The differential reactivity worths of the Mo sleeve and stainless -steel clad 
cadmium (Cd) sleeve as obtained from one -dimensional S n calculations were not in 
agreement with measured reactivity worths. 

In summary, then, Mo cross sections from the ENDF/B #1111 evaluation are ade- 
quate for predicting the absolute multiplication factor of the seven reactors presented in 
this report. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 4, 1972, 

112-02. 
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